Spin-state transition in LaCoOa: direct neutron spectroscopic 
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Abstract 

A gradual spin-state transition occurs in LaCoOs around T ~ 80 — 120 K, whose detailed nature 
remains controversial. We studied this transition by means of inelastic neutron scattering (INS), 
and found that with increasing temperature an excitation at ~ 0.6 meV appears, whose intensity 
increases with temperature, following the bulk magnetization. Within a model including crystal 
field interaction and spin-orbit coupling we interpret this excitation as originating from a transition 
between thermally excited states located about 120 K above the ground state. We further discuss 
the nature of the magnetic excited state in terms of intermediate-spin (IS, t^e^, S = 1) vs. high- 
spin (HS, t| fl e 2 , S = 2) states. Since the ^-factor obtained from the field dependence of the INS is 
g ~ 3, the second interpretation looks more plausible. 
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Due to its rich and in many respects puzzling properties, LaCoC>3 keeps attracting atten- 
tion and remains a controversial topic. It is known that the ground state is nonmagnetic, 
corresponding to a low-spin (LS) state of Co 3+ ions (t| g , 5 = 0). However, with increasing 
temperature (as well as with La— >Sr substitution) first a crossover into a magnetic, but still 
insulating state appears at about 80-120 K, followed by another crossover into a "bad metal- 
lic" , magnetic state at T ~ 400 — 600 K. The original interpretation of the low-temperature 
crossover was done in terms of thermally-induced population of the low-lying high-spin (HS) 
state Q; this process is furthermore favorized by thermal expansion, since the HS Co« has 
much larger radius (~0.75 A) than the LS state (~0.685A). Later, especially after LDA+U 
band structure calculation have become available Q], another interpretation was put for- 
ward: within this scenario, the first crossover at ~ 100 K would be due to a transition into 



an intermediate-spin 



mental evidences 



S) state. This interpretation was supported by a number of experi- 



Q,y,y,H 



111 ]. None of these arguments however gave a definite 
proof that the first thermally-excited state is indeed the IS one. Very recent measurements 



e. very recent mea; 



indicate that the first excited state could still be the HS state 

Also theoretically the situation is not clear: Hartree-Fock calculations showed that the HS 
state, or the HS-LS ordered state is more stable than the IS state [17], in contrast to LDA+U 
calculations |2j. Thereby, model calculations on a C0O6 cluster explicitly including the Co- 
O hybridization can not reproduce an IS ground-state [18]], indicating that the proposed 
mechanism why LDA+U finds an IS as first excited state, namely large covalency, is rather 
questionable. 

With this controversy in mind, we undertake a neutron scattering study of LaCoC>3 at 
different temperatures with the goal of identifying the energy level of the thermally excited 
state of Co 3+ . Indeed, we discovered that a rather unusual feature in the spectrum appears 
with increasing temperature in forms of thermally-induced relatively sharp inelastic peak 
at an energy-transfer of ~ 0.6 meV Q]. The intensity of this peak strongly increases with 
T following the behavior of the magnetic susceptibility xCO> suggesting that the inelastic 
scattering occurs between thermally populated magnetic states of LaCoC>3. The position and 
the temperature dependence of the intensity of this peak also coincide with the excitations 
observed in LaCoOs by ESR 2, l^of - By analyzing the features of this novel excitation, 
and combining it with model calculations, we discuss the two possible scenarios mentioned 
above. 
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FIG. 1: (Color online) Temperature evolution of the INS profiles measured in LaCo03. The filled 
circles correspond to the LaA103 nonmagnetic reference compound at T = 50 K. The lines are the 
result of least-squares fits using Gaussian functions to describe the lineshape of the transition. For 
clarity, an offset has been added to the various curves. 

Polycrystalline LaCoOa was prepared by standard sintered techniques using La203 and 
C03O4 of a minimum purity of 99.99%. The respective amounts of starting reagents were 
mixed and calcinated at temperatures 1000-1200° C during at least 100 h in air, with several 
intermediate grindings. The sample was checked by x-ray diffraction and found to be single 



phase within experimental accuracy. The space group R3c and lattice parameters o 
5.4433(1) A, c = 13.0932(4) A are in agreement with previously published data 
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The inelastic neutron scattering (INS) measurements were performed on the time-of-flight 
spectrometer FOCUS 



23| installed at the spallation neutron source SINQ at Paul Scherrer 
Institut, Villigen, Switzerland. Zero field experiments were carried out in the temperature 
interval of 1.5 — 100 K using a conventional helium cryostat. The data were collected using 
incoming neutron energies of 3.5 and 20 meV, giving an energy resolution at the elastic 
position [full width at half maximum (FWHM)] of 0.1 and 1.6 meV, respectively. The triple- 
axis spectrometer TASP with final neutron energy 4.7 meV was used for the measurements 
in external magnetic field up to H = 6 T. 

The high-resolution low-energy transfer inelastic spectra for a few selected temperatures 
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FIG. 2: Observed (circles) and fitted (solid line) temperature dependence of the integral intensity 
of the INS peak at 0.6 meV. The resulting level scheme is shown in the inset. 

are shown in Fig. ^ There are no excitations in the energy window E < 1.5 meV for 
temperatures T < 30 K. A single inelastic peak at an energy transfer SE = 0.61 ±0.05 meV 
was found at intermediate temperatures starting from T ~ 30 K. A strong broadening of the 
transition was observed with increasing temperature. Note that the spectra obtained from 
a non- magnetic reference compound, LaAlOa, remain structureless at all temperatures. The 
high-energy transfer spectra observed for LaCoOs exhibit several broad inelastic peaks at 
about 10, 14 and 22 mev (not shown). However, all these peaks exhibit clear increase of 
their intensity with scattering vector and temperature. Therefore we conclude that they are 
due to phonon scattering, in agreement with previously published data No evidence for 
other magnetic excitations was observed in the LaCo03 spectra within the available energy 
window. 

For noninteracting ions the thermal neutron cross section for the transition \Ti) —>■ \Tj) 



is given in the dipole approximation by 



24] 



1 



exp 



Ei 



F 2 (Q)\(T j \J ± \T i )\ 2 5(E i - Ej ± hu) 



(1) 



dVtduj Z 1 V kftT 
Here F 2 (Q) is the magnetic form factor, J± is the component of the total angular momentum 
operator perpendicular to the scattering vector Q, and Z is the partition function. It follows 



from Eq. (0) that the energy gap AE = Ei — E between the ground-state and the excited- 
state can be deduced either directly from the position of the corresponding inelastic peak 
(in case of nonzero matrix element |(r | Jj_|rj)|), or from the temperature dependence of 
the transition between two excited levels 8E = Ej — Ei which is governed by Boltzman 
statistics. Note that the direct transition AE out of the ground-state was observed neither 
in the previous INS expend QQ, nor in our current n—ents, ,nost likely due to 
selection rules. Therefore, in order to determine the energy of the excited state we apply the 
least-squares fitting procedure to the temperature dependence of the integrated intensity I 
of the INS signal as shown in Fig. 121 The position of the excited states turns out to be 



10.3 ± 1 meV, which coincides well with the results obtained from ESR (12 meV, ref. 



Our estimation is based on a temperature independent level splitting scheme. Although 
we can not exclude a slight variation of the position of the excited states due to thermal 
expansion of the unit cell, a level crossing of the ground- and excited states as suggested from 
the LDA+U calculations Q can be excluded, since this would result in a non-monotonic 
temperature behavior of / around T < 80 K, which was not observed in our experiment. 
Furthermore, the position of the peak at 5E = 0.6 meV is unaltered, suggesting that the 
trigonal CF remains nearly constant in this temperature range. Thus, our results imply 
that the first broad peak in magnetic susceptibility at ~ 100 K is due to a gradual thermal 
population of the excited levels rather than a modification of the level scheme due to a phase 
transition. 

The observed magnetic INS, which was obtained as the difference of the intensities at 50 
and 5 K, is shown in Fig. El A clear shift of the transition to the higher energy ~ 1.5 meV 
was observed in magnetic field H = 6 T compared to the zero-field spectrum, thus firmly 
establishing its magnetic origin. Due to weakness of the signal we can not conclude whether 
the peak is split in the external magnetic field. The change in energy of this peak from 
0.6 meV to about 1.5 meV in magnetic field of 6 T is in good agrement with the (/-factor 
measured from ESR experiments 12J. 

Let us discuss the possible origin of this excited state. There are two possibilities: either 
high-spin or a intermediate-spin states of Co 3+ . The HS state with 5* = 2 has, in a cubic 
CF, the occupation t\ g e 2 g1 i.e., it has half-filled shell of t\ g e 2 g1 say with spins up, and one extra 
spin-down electron on a triple-degenerate t2„-level, which can be described by an effective 

n 

orbital moment L = 1 [26]. Total multiplicity of this state is (2S + 1)(2L + 1) = 15. Spin- 
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FIG. 3: (Color online) The magnetic inelastic scattering at T = 50 K in T (circles) and 6 T 
(boxes) applied field. 

orbit coupling splits this state into the lowest-lying triplet J = S — L — 1, next is a quintet 
J = 2, and the highest-lying is state has J = 3. If the system is strongly distorted there 
will be a ground-state orbital singlet and a higher excited orbital doublet. In the left panel 
of Fig. 0]we show the energy level diagram for the high-spin state as a function of trigonal 
distortion. This energy level diagram has been calculated for a CoOg - cluster, including 
full multiplet theory, spin-orbit coupling and Co-0 hybridization. For the Slater integrals 
and the spin orbit coupling atomic Hartree-Fock values are used, the hopping parameters 
are according to Harrisons's rules j^. The calculations have been done with the use of the 



program XTLS8 



281 ] . For the HS there are two places in the energy- level diagram, where 



an excitation of 0.6 meV can happen. If the trigonal crystal field is relatively small, the 
J = 1 triplet will be split by this crystal field. On the other-hand, if it is rather large, the 
orbital singlet with S = 2 will be split due to spin-orbit coupling. In both cases the splitting 
is governed by second order effects and the resulting splitting is much smaller than the 
perturbing interaction. The scenario of a small crystal-field with respect to the spin-orbit 
coupling has been discussed in quite details recently Q]. The scenario of a large crystal- field 
with respect to the spin-orbit coupling is equivalent to a spin only scenario. 

One should also consider what would be the situation if the first excited state is an 
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FIG. 4: Schematic diagram of the excitation spectrum of high-spin (left) and intermediate-spin 
(right) Co 3+ as a function of trigonal distortion (see text). 



IS Co 3+ , which follows from LDA+U theoretical calculations M and which was used to 
interpret a number of experimental data llll ]. This case is actually 

much more interesting and more difficult to treat theoretically. First thing to note is that 
in this case we have one electron in the e 9 -shell and one hole (t\ g occupation) in the t2 9 - 
shell. Their Coulomb interaction strongly depends on the particular orbital occupation. 
The (x 2 — y 2 )-electron has strong attraction to the (xy)-hole, so that (x 2 — y 2 ) 1 {xy) 1 state 
has much lower energy than e.g. (z 2 ) l (xy) 1 . Thus we can consider the lowest states on 
a basis of (x 2 — y 2 ) 1 (xy) 1 , (x 2 — z 2 ) l (xz) 1 and [xy 2 — z^i^yz) 1 [3]. Therefore the total 
orbital degeneracy of the IS state in a cubic CF is 3 and not 6 {3t2 g x 2e g ) as one could have 
expected. We can thus again describe these states by the effective orbital triplet L = 1; but 
because of the more complicated type of the basis states, the maximum magnitude of the 
magnetic moment is not 1 but \ (i.e. L z = {|,0, — \} 29].). In other words, one has an 
effective orbital (/-factor of |. As a result, we are dealing with 9 states (2L + l)(25 l + 1) 
(with S = 1 for IS state), which are split by the spin-orbit coupling into multiplets with 
J = 2, 1, 0. However, in this case the quintet J = 2 is the lowest state. Thus, the multiplicity 
of the IS state in cubic CF is 5. Strong enough distortions, or orbital ordering, modify the 
energy-level scheme as shown in the right panel of Fig. HJ If the distortions are larger than 
the spin-orbit coupling constant, the ground-state becomes an orbital singlet. This orbital 
singlet is split due to second-order spin-orbit interactions into two levels that could very well 
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be 0.6 meV apart from each other. In cubic symmetry, the J = 2 quintet originating from 
the IS state is however also split due to second order interactions. This splitting seems to 
be somewhat larger than the measured value of 0.6 meV. 

If we now compare our results with other measurements, we notice that from a comparison 
of the magnetic susceptibility with the anomalous expansion coefficient of LaCo03, Zobel et 
al. [n| concluded that the degeneracy of the first magnetic excited state is 3. This leaves 
only two possibilities open. The first magnetic excited state in LaCo03 can be a HS state 
with a small non-cubic crystal field, or it can be an IS state with a large non-cubic crystal 
field or orbital ordering. There is one striking difference between these scenarios: this is the 
predicted q-factor. The HS state with a small distortion is a triplet with a g-factor of about 
3.5 3, 13 whereas the IS with strong distortion is a triplet with a ^-factor of about 2.0. 
ESR measurements found a ^-factor of 3.35 - 3.55 12j ]. supporting the HS state; this also 
agrees with our results since we obtained g ~ 3 (see Fi g. El) . On the other hand a fit to the 
magnetic susceptibility yields a ^-factor of about 2.28 supporting the IS scenario. 

To summarize, we observe a novel inelastic excitation in LaCo03 which is due a thermally 
excited magnetic state of Co 3+ ions. This confirms the presence of thermally induced spin- 
state transition (or rather crossover) at T ~ 100 K from the LS Co 3+ to a magnetic HS or IS 
state. We discuss both possibilities theoretically and show that one can explain both our and 
other results (thermodynamic, ESR) in the framework of a HS-triplet as first excited state, 
with the g-factor ~ 3.5, weakly split by small distortions from the cubic symmetry. Another 
possibility would be the IS state with orbital ordering or strong non-cubic crystal-fields, 
which however would result in a spin-only system with a g-f actor ~ 2.0, a value difficult to 
reconcile with our experimental data. Thus the first interpretation (HS excited state) seems 
to us more plausible. 
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